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Abstract 

Inclusive non-diffractive photoproduction of p(770)°, K*(892)° and 0(1020) mesons is 
investigated with the HI detector in ep collisions at HERA. The corresponding average jp 
centre-of-mass energy is 210 GeV. The mesons are measured in the transverse momentum 
range 0.5 < pr < 7 GeV and the rapidity range \yi a b\ < 1- Differential cross sections 
are presented as a function of transverse momentum and rapidity, and are compared to the 
predictions of hadroproduction models. 
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1 Introduction 



High energy particle collisions, which give rise to large multiplicities of produced hadrons, pro- 
vide an opportunity to study the hadronisation process, whereby the quarks and gluons produced 
in the initial interaction become colourless hadrons. Since most of these hadrons are produced 
at low values of transverse momentum, perturbative quantum chromodynamics (pQCD) is not 
applicable to this process, which is described instead using phenomenological models, the most 
successful of which are the string [1] and the cluster [2] fragmentation models. These can pro- 
vide a reasonable description of the hadronisation process provided the many free parameters 
they contain are tuned to the data. 

The production of long-lived hadrons and resonances at high energies has been studied in 
detail in electron-positron (e + e~) collisions at LEP using Z° decays [3]. Measurements in high 
energy hadronic interactions have so far been restricted to long-lived hadrons and hadrons con- 
taining heavy quarks. Recently, the production of the hadronic resonances p(770)°, i^*(892)° 
and 0(1020) has been measured in heavy-ion and proton-proton (pp) collisions at RHIC [4]. 
The electron-proton (ep) collider HERA allows the study of particle production in quasi-real 
photon-proton (7^) collisions. The comparison of RHIC and HERA results is of particular in- 
terest, since the nuclear density at HERA is much lower than that at RHIC while the 7p and 
nucleon-nucleon collision energies are similar. 

In this paper, measurements of the inclusive non-diffractive photoproduction of the reso- 
nances p(770)°, if* (892)° and 0(1020) at HERA are presented for the first time. The measure- 
ments are based on the data recorded with the HI detector during the year 2000, when positrons 
of energy 27.6 GeV collided with 920 GeV protons at an ep centre-of-mass energy of 319 GeV, 
providing on average a 7p centre-of-mass energy of (W) = 210 GeV. The data correspond to 
an integrated luminosity of C = 36.5 pb" 1 . 



2 Phenomenology and Monte Carlo Simulation 

The HI coordinate system has as its origin the position of the nominal interaction vertex. The 
outgoing proton beam direction defines the positive z-axis and is also referred to as the "for- 
ward" direction. The polar angle 9 is defined with respect to this direction. The pseudorapidity 
is given by r]i ab = — ln(tan(6 l /2)). The laboratory frame rapidity yi ab of a particle with energy 
E and longitudinal momentum p z is given by yi ab = 0.5 \n[(E + p z ) /(E — p z )]. 

The invariant differential cross section for meson production can be expressed as a function 
of the meson's transverse momentum p T and its rapidity y iab , assuming azimuthal symmetry. 
Hadrons produced in hadronic collisions are approximately uniformly distributed in the central 
rapidity range, while their transverse momentum spectra fall steeply with increasing p?- It is 
convenient to parametrise the invariant differential cross section of the produced hadrons with 
a power law distribution, 

1 d 2 a^ p _ A 

n dp 2 T dy lab = {E To + E* m y ' (1) 
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where EJp n = \J + p\ — m Q is the transverse kinetic energy, m is the nominal reso- 
nance mass, A is a normalisation factor independent of px and E To a free parameter. When 
E!p n < E To , the power law function (1) behaves like a Boltzmann distribution exp(— E^ m /T), 
with T = Ex /n. This exponential behaviour of hadronic spectra follows from a thermody- 
namic model of hadroproduction [5]. In this framework, the parameter T plays the role of the 
temperature at which hadronisation takes place. At high E!p n , the power law originates from a 
convolution of the parton densities of the colliding particles with the cross sections of parton- 
parton interactions. The normalisation coefficient A is related to the single differential cross 
section da/dyi ab obtained after the integrating equation (1) over p\: 

A= da_ (n-l)(n-2)(E To )"- 1 _ 
dy iab 2tt(E To + (n- 2)m ) 

Monte Carlo calculations are used both to correct the data and in comparisons with the measure- 
ments. Direct and resolved photoproduction events are simulated using the PYTHIA [6] and 
the PHOJET [7] Monte Carlo generators. In both cases, the hadronisation is based on the string 
fragmentation model [8]. For data corrections, the parameter settings obtained by the ALEPH 
collaboration [9] are used for the fragmentation of partons. The effects of Bose-Einstein corre- 
lations (BEC) on the invariant mass spectra of like-sign and unlike-sign pion pairs are included 
using a Gaussian parametrisation of the correlation function [9]. The photoproduction events 
generated using PYTHIA and PHOJET are passed through the simulation of the HI detector 
based on GEANT [10] and through the same reconstruction and analysis chain as used for the 
data. 



3 Experimental Conditions 
3.1 HI Detector 

The HI detector is described in detail elsewhere [11]. A brief account of the components that 
are most relevant to the present analysis is given here. 

The ep interaction region is surrounded by two large concentric drift chambers (CJCs), 
operated inside a 1.16 T solenoidal magnetic field. Charged particles are measured in the pseu- 
dorapidity range —1.5 < r]i ab < 1.5 with a transverse momentum resolution of o PT jp T « 
0.005 ■ p T /GeV © 0.015 [12]. The specific energy loss dE/dx of the charged particles is 
measured in this detector with a relative resolution of 7.5% for a minimum ionising track [13]. 

A finely segmented electromagnetic and hadronic liquid argon calorimeter (LAr) covers the 
range— 1.5 < rji ab < 3.4. The energy resolution of this calorimeter is a(E)/E = 0.11/ y^E/G&V 
for electromagnetic showers and a(E)/E = 0.50/ \/E/GeV for hadrons as measured in test 
beams [14]. 

Photoproduction events are selected with a crystal Cerenkov calorimeter (positron tagger) 
located close to the beam pipe at z = —33.4 m, which measures the energy deposited by 
positrons scattered at angles of less than 5 mrad. Another Cerenkov calorimeter, located at 
z = —103 m (photon tagger), is used to determine the luminosity by measuring the rate of 
photons emitted in the Bethe-Heitler process ep — > epr). 
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3.2 Event Selection 



Photoproduction events are selected by a trigger which requires a scattered positron to be mea- 
sured in the positron tagger, an event vertex determined from charged tracks and three or more 
charged tracks reconstructed in the CJCs, each with transverse momentum px > 0.4 GeV. The 
photon virtuality Q 2 is smaller than 0.01 GeV 2 , due to the positron tagger acceptance. The pho- 
ton energy is determined from the difference between the positron beam energy and the energy 
measured in the positron tagger. 

In order to reduce the non-ep background and to ensure good reconstruction of the event 
kinematics, the following criteria are applied: 

• Events are selected if the reconstructed ^yp centre-of-mass energy lies within the interval 
174 < W < 256 GeV for which good positron detection efficiency is established. This 
corresponds to an average ^p centre-of-mass energy of {W) = 210 GeV. 

• Events are rejected if a photon with energy E 1 > 2 GeV is detected in the photon tag- 
ger. This suppresses the background arising from random coincidences of Bethe-Heitler 
events in the positron tagger with beam-gas interactions in the main HI detector. 

• Events are selected if the z coordinate of the event vertex, reconstructed using the CJCs, 
lies within 35 cm of the mean position for ep interactions. 

Background from elastic and diffractive events is suppressed by the above trigger require- 
ments. To further reduce the contribution of diffractive processes, the presence of an en- 
ergy deposit of at least 500 MeV is required in the forward region of the LAr, defined by 
2.03 < r]iab < 3.26. Monte Carlo studies show that, with this requirement, less than 1% of the 
final event sample consists of diffractive events with Xp < 0.05, where Xp = M\/W 2 and 
Mx is the invariant mass of the diffractive system. 

In total, about 1.8 x 10 6 events satisfy the above selection criteria. 

3.3 Selection of p(770)°, K*(892)° and 0(1020) Mesons 

The mesons are identified using the p(770)° -> tt+tt^, (892)° -> K + ir~ or iT (892)° -> 
K~7r + and 0(1020) — > K + K~ decays 1 . Charged tracks reconstructed in the CJCs withp^ > 
0.15 GeV and pseudorapidity \r][ a b\ < 1.5 are considered as charged pion or kaon candidates. 
Since most of the charged particles in ep collisions are pions, no attempt to identify pions is 
made, while identification criteria for charged kaons are applied for the extraction of the K*° 
and (ft signals. This is done by measuring the momentum-dependent specific energy loss dE/dx 
in the CJCs. This method gives a significant improvement in the signal-to-background ratio 
for low pr mesons, pt < 1.5 GeV, where the dE/dx resolution allows good particle identi- 
fication. For high mesons, p? > 1.5 GeV, the dE/dx method is inefficient and therefore 

'in the following, the notation K*° is used to refer to both the K*° and K*° mesons unless explicitly stated 
otherwise. 
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particle identification criteria are not applied. Such tracks are considered as both pion and kaon 
candidates and their four-momenta are determined from the track measurements using the cor- 
responding mass hypothesis [15]. Vector meson candidates are reconstructed from these four- 
momenta. The kinematic range for the reconstructed neutral mesons is restricted to \yi a b\ < 1 
andp T > 0.5 GeV. 

To extract the p°, K*° and 4> signals, the distributions of respective invariant masses of their 
decay products, m^+^-, m K ± 7T ^ and m K + K ~, are fitted using a function composed of three 
parts: 

F(m) = B(m) + i?(m) + ^S(m). (3) 

The terms correspond to contributions from the combinatorial background, B(m), from reflec- 
tions which result from decays other than the signal under consideration, R(m), and from the 
relevant signal, S(m), respectively. 

The combinatorial background function is taken to be: 

B{m) = (do + aim + (12m 2 + a 3 m 3 ) ■ B°(m) , 

where a , a 1? a 2 and a 3 are free parameters, and B°(m) is the invariant mass distribution of the 
like-sign charged particle combinations: iv ± 7r ± for the p° and K ± ir ± for the A'* . The shape of 
the combinatorial background for <p is described by the following function: 

B(m) = him 2 - 4m 2 K ) b2 e- bsm , 

where b\, b 2 and 63 are free parameters and mx is the kaon mass. 

The second term, R(m), in (3) represents the sum of the reflections; for example, charged 
particles from the decay A"* — ► K ± tt t with the kaon misidentified as a charged pion will 
give rise to structure in the m^ n - spectrum and must be taken into account as a separate 
contribution. In addition, there are two other contributions to the m^^- spectrum in the mass 
region of interest. These arise from the decays a>(782) — > n + ix~ and lj(782) — > iv + n~7r in 
which the ir° is not observed. For the c<j(782) meson, the production rate relative to that of the 
p° is varied within the range 1.0 ± 0.2, which is consistent with measurements of the cj(782) / p° 
ratio in hadronic collisions [16] and in Z° boson decays [17]. The c<j(782) branching ratios 
are taken from [15]. The five major reflections in the m K ±^ spectrum are due to: the decay 
p° — > n + ix~ with the n + or n~ misidentified as a charged kaon; the decays c<j(782) — > ix + ix~ and 
u;(782) — > 7r + 7r~7r° with the 7r° not observed and with one of the n + or n~ mesons misidentified 
as a charged kaon; the decay <\> — > K + K~ with one of the kaons misidentified as a charged 
pion and a self-reflection from the K*°, where the pion and kaon are interchanged. For the 
m K+K- spectrum, there are no reflections from known resonances in the invariant mass region 
of interest. Therefore, the shapes of the reflections are taken from Monte Carlo calculations. 
The contribution of the reflections from the p°, K*° and mesons is tied to the production rates 
determined in this analysis and is therefore iteratively calculated. 

The function S(m) used to describe the signal in (3) is a convolution of a relativistic Breit- 
Wigner function BW(m) and a detector resolution function r(m, m'). The relativistic Breit- 
Wigner function 

BW( A m mT(m) 

(m 2 — ml) 2 + ml T 2 {m) ' 



is used with 




) 



2Z+1 



r(m) = r 



m 



where A is a normalisation factor, T is the resonance width, / = 1 for vector mesons, m is the 
resonance mass, q is the momentum of the decay products in the rest frame of the parent meson, 
and g is their momentum at m = m . The cross sections cited in this paper assume that the 
meson signal is defined as the integral of the relativistic Breit-Wigner function (4) in the region 
±2.5r around the mass m . Monte Carlo studies show that a non-relativistic Breit-Wigner 
function with width r res provides a good description of the detector resolution function: 



For the K analysis, the resolution parameter is determined from Monte Carlo with r res = 
12 MeV. It is small compared to the width of the K*° meson (50.3±0.6 MeV) [15], leading only 
to a small change in the shape of the resonance. For the 0, T res is comparable to the width of 
the meson (r = 4.26 ± 0.05 MeV) [15]. As a result, the shape of the <p signal is significantly 
changed, and hence the detector resolution r res is taken as a free parameter in the fit. It is found 
to vary from 3.4 MeV to 6.0 MeV, increasing with the px of the (ft meson. 

For the p° meson, the detector resolution is significantly smaller than its width. However, 
BEC between the p° decay pions and other pions in the event strongly distort the p° line shape. 
The BEC plays an important role in broadening the p° mass peak and in shifting it towards 
lower masses. Similar effects are observed in pp and heavy-ion collisions at RHIC [4] and in 
e + e~ collisions at LEP using Z° decays [17]. It is therefore important to check that the Monte 
Carlo model used for the extraction of the cross sections describes the di-pion spectra in the 
data. The data spectra and the Monte Carlo simulations with and without BEC are shown in 
figure 1 . The Monte Carlo model with BEC is in a good agreement with the data in the region of 
the p° resonance, whereas the model without BEC fails to describe the di-pion mass spectrum. 

The results of fitting the function (3) to the m n + n - data in the mass range from 0.45 to 
1.7 GeV with the contributions due to the combinatorial background and the reflections are 
shown in figure 2a), and after combinatorial background subtraction in figure 2b). In this mass 
range, the signal from the Kg, /o(980) and / 2 (1270) mesons is taken into account. The Kg sig- 
nal is fitted using a Gaussian centred on the nominal mass and with fixed width. The relativistic 
Breit-Wigner function given in equation (4) is used for the / (980) and / 2 (1270) mesons. In 
the fit, the resonance masses m and the yields are free parameters. The p° and / 2 (1270) widths 
are fixed to the Particle Data Group [15] values and the / (980) width is fixed to 70 MeV. Due 
to the small signal and the non-trivial background behaviour, which lead to large uncertainties, 
cross sections for / (980) and / 2 (1270) meson production are not measured here. 

The A'* signal is measured under the assumption that there is no difference between the 
particle and antiparticle production rates, and the signal obtained from the m K ± 7r ^ spectrum is 
divided by 2 to determine the K*° rate in the following. The result of fitting the function (3) 
to the m K ± 7r ^ data in the mass range from 0.7 to 1.2 GeV with the contributions due to the 
combinatorial background and the reflections is shown in figure 2c). In the fit, the K*° width is 
fixed to the nominal value while the mass parameter is left free. The result for the K*° mass is 
compatible with the world average [15]. 



r(m, mf) 



re* 



(5) 



2tt (m - m') 2 + (T res /2) 2 ' 
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The result of fitting function (3) to the m K + K - data in the mass range from 0.99 to 1.06 GeV, 
together with the background contribution, is shown in figure 2d). In the fit, the width, T , is 
fixed to the nominal value while the mass is left a free parameter and is found to be compatible 
with the world average [15]. 

3.4 Cross Section Determination and Systematic Errors 

The invariant differential cross section for p°, A'* and <f) meson production is measured in 
the yi a b region from —1 to 1 in seven bins in transverse momentum from 0.5 to 7 GeV. It is 
calculated according to: 

1 cPa-v _ N 

n dp% dy iab 7r ■ C ■ BR ■ <3> 7 ■ e • Ap\ ■ Ay lab ' 

where N is the number of mesons from the fit in each bin. The corresponding bin widths are 
Ay iab = 2 and Ap^ = 2p^ n ApT- Bin centre corrections based on equation (1) are applied 
to define the value of p^ n at which the differential cross section is measured. C denotes the 
integrated luminosity and e the efficiency. The branching fractions BR are taken from [15] and 
are equal to 1, 0.67 and 0.49 for p° — > ix + n~ , K*° — > K^ir^ 1 and </> — > K + K~ , respectively. 
The photon flux $ 7 = 0.0127 is calculated using the Weizsacker- Williams approximation [18]. 

The single differential cross section for p°, K*° and meson production for pt > 0.5 GeV 
is measured in four bins in rapidity from —1 to 1 according to: 

da^ p _ N 

dyi a b £ ■ BR ■ $ 7 • e • Ay iab ' 

Here, the bin width is Ayi a b = 0.5. 

The fit procedure described in the previous section is repeated to determine the number of 
mesons, N, in each measurement bin, calculated as an integral over the signal function (4) 
within ±2.5r around the meson mass. Similarly, the total visible cross section for p°, K*° and 
(j) meson production is measured from the number of mesons fitted in the range \yi a b\ < 1 and 
p T > 0.5 GeV. 

The efficiency is given by e = e rec ■ A eta g • ^3 ■ ttrig- The reconstruction efficiency for the 
mesons, e rec , includes the geometric acceptance and the efficiency for track reconstruction. It 
is calculated using Monte Carlo data and is at least 45% at low pt and rises to about 90% with 
increasing p?- For the acceptance determination, the Monte Carlo generators are reweighted 
to model the observed ^-dependences . The average acceptance of the positron tagger, A e t ag , 
is about 50%, as determined in [19]. The trigger acceptance, A3, arises from the requirement 
that at least three tracks are reconstructed in the CJCs with p T > 0.4 GeV. It is determined 
from Monte Carlo simulations with PYTHIA and PHOJET and varies from 50% to 95%. The 
trigger efficiency, e tr i g , is calculated from the data using monitor triggers. It is about 90%. The 
efficiencies and acceptances as calculated from the PYTHIA and PHOJET simulation are found 
to be consistent. Small residual differences, attributed to different track multiplicity predictions, 
are taken into account in the systematic uncertainties of the measurement. 
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The statistical error varies from 7 to 15% for the p°, 10 to 18% for the K*° and 13 to 24% 
for the meson cross sections. The systematic errors arise from the uncertainties in the track 
reconstruction efficiency (4%) and the trigger efficiency (up to 6%), the variation of the / (980) 
width by ±30 MeV in the p° fit (up to 7%), the uncertainties in the dE/dx kaon identification 
procedure (6% for the K*° and 12% for the 0) and the luminosity calculation (2%), the variation 
of the background shape (5%) and the variation of the assumptions about the normalisation of 
the contributions from the reflections (4% for the p° and up to 15% for the K*°). The total 
systematic error varies from 10 to 12% for the p°, 11 to 21% for the K*° and 10 to 17% for the 
meson cross sections. 



4 Results and Discussion 

The inclusive non-diffractive photoproduction cross sections for p°(770), K*°(892) and 0(1020) 
mesons in the kinematic region Q 2 < 0.01 GeV 2 , 174 < W < 256 GeV, and for p T > 0.5 GeV 
and \yiab\ < 1 are found to be: 

al P ls (lP P° x ) = 25600 ± 1800 ± 2700 nb; 

^L(lP~^ K*°X) = 6260 ± 350 ± 860 nb; 

aJf s (7p-*0X) = 2400 ± 180 ± 340 nb. 

The first error is statistical and the second systematic. Note that the K*° cross section is the 
sum of the particle and antiparticle contributions divided by 2. 

The differential cross sections for the photoproduction of p°, K*°, and mesons are pre- 
sented in tables 1 and 2 and in figure 3. Within the rapidity range of this measurement, the 
resonance production rates are constant as a function of rapidity, within errors. The transverse 
momentum spectra of the p°, K*° and mesons can be parametrised by function (1), where 
da/dyi a b in equation (2) corresponds to the average value of the cross section over central ra- 
pidities, (dcr/dyiab)^^^. In the fit, the value of the power n is fixed to be 6.7, as derived 
previously from measurements of charged particle spectra by the HI collaboration [20] which 
gave n = 6.7 ± 0.3. The power law distribution, with this value of n, describes K s meson, A 
baryon [21] and D*^ 1 meson production [22] at HERA, as is shown in figure 4. A similar shape 
of the transverse momentum distribution, but with different values of the parameters n and E To , 
was reported for charged particles produced in hadronic collisions [23]. The results of the fits of 
the data to function (1) are shown in figure 3a). In table 3, the parameters of the fit and the av- 
erage transverse kinetic energy (E!p n ), the average transverse energy (E T ) = (E!p n ) + m and 
the average transverse momentum (p T ) = a/ (E t ) 2 — ml derived from (1) are presented. The 
errors include the experimental uncertainty on the value of n. Also given are the (pr) values 
measured at RHIC in pp and Au-Au collisions [4] . 

It is interesting to observe that the resonances with different masses, lifetimes and strangeness 
content are produced with about the same value of the average transverse kinetic energy (E^ m ). 
This observation supports the thermodynamic picture of hadronic interactions [5], in which 
the primary hadrons are thermalised during the interaction. The values of (p T ) for p°, K*° 
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and mesons are similar in ^p and pp collisions with about the same centre-of-mass energy 
y/s ~ 200 GeV, while these values are all higher in Au-Au collisions. 

The PYTHIA and PHOJET models do not describe the shape of the measured pt spectra. 
Moreover, contrary to the data, the Monte Carlo px spectra are not described by the power law 
function (1). These observations are illustrated in figures 3c) and 3d). 

The measurements in the visible kinematic range of the p°, K*° and <p mesons, pr > 
0.5 GeV and \yi ab \ < 1, are extrapolated to the full px range using the parametrisation (1) 
to determine the total inclusive non-diffractive photoproduction cross sections. The extrapola- 
tion factors are of order two. In the rapidity interval \yi a b\ < 1 an d integrated over the full pr 
range the following cross section ratios are obtained: 

R(K*°/p°) = 0.221 ±0.036; 

R(<P/p°) = 0.078 ±0.013; 

R(<p/K*°) = 0.354 ±0.060. 

The errors are given by the statistical and systematic errors added in quadrature. PYTHIA and 
PHOJET, with the strangeness suppression factor \ s = 0.286 [9], predict the ratios 0.200, 0.055 
and 0.277, respectively, which are similar to the measured values, but are all somewhat lower 
than these. 

In table 4, R((j)/K*°) is compared to the corresponding ratios measured by STAR in pp and 
Au-Au collisions [4] at ^snn = 200 GeV. Although the rapidity ranges at the HI and RHIC 
experiments differ 2 , the resulting ratios for pp and ^p interactions are very close. However, the 
corresponding result in Au-Au collisions is observed to be higher. 

5 Conclusions 

First measurements of the inclusive non-diffractive photoproduction of p(770)°, K* (892)° and 
0(1020) mesons at HERA are presented. The differential cross sections for the production of 
these resonances as a function of transverse momentum are described by a power law distri- 
bution while the single differential cross sections as a function of rapidity are observed to be 
flat in the visible range. Despite their different masses, lifetimes and strangeness content, these 
resonances are produced with about the same value of the average transverse kinetic energy. 
This observation supports a thermodynamic picture of hadronic interactions. 

The description of the shape of the p° resonance produced in 7p collisions at HERA is 
improved by taking Bose-Einstein correlations into account. A similar effect is observed in pp 
and heavy-ion collisions at RHIC and in e + e~ annihilation at LEP, using Z° decays. 

The cross section ratios R(K*°/p°), R(4>/p°) and R((p/K*°) are determined, and R(<p/K*°) 
is compared to results obtained in pp and heavy-ion collisions by the STAR experiment at RHIC. 
The ratio R((p/K* ) measured in 7p interactions is in agreement with the pp results, while this 
ratio is observed to be smaller than the result obtained in Au-Au collisions. 

2 The difference in rapidity between the laboratory frame and the jp frame is about two units at HI. 
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2 

- 3 2? [nb/(GeV) 2 ] 

vr dp^dy lab 


Pt iyj£y J 


bin 

Pt 


P 


{IS. -+- J\ ) / z 


9 


[0.5,0.75] 


0.63 


5610 ± 870 ± 590 


1190 ± 130 ± 200 


383 ± 54 ± 60 


[0.75,1.0] 


0.87 


2440 ± 180 ± 260 


621 ± 68 ± 80 


264 ± 34 ± 37 


[1.0,1.5] 


1.22 


680 ± 55 ± 70 


176 ± 18 ± 21 


76 ± 12 ± 11 


[1.5,2.0] 


1.72 


142 ± 15 ± 15 


48.0 ± 5.2 ± 5.1 


19.1 ± 3.3 ± 1.9 


[2.0,3.0] 


2.41 


29.9 ± 2.3 ± 3.1 


8.96 ± 0.90 ± 0.98 


3.48 ± 0.76 ± 0.34 


[3.0,4.0] 


3.43 


3.06 ± 0.42 ± 0.33 


1.21 ± 0.17 ± 0.14 


0.46 ± 0.11 ± 0.08 


[4.0, 7.0] 


5.09 


0.276 ±0.037 ±0.033 


0.079 ±0.014 ±0.009 


0.0335 ±0.0081 ±0.0057 



Table 1: Inclusive non-diffractive photoproduction invariant differential cross sections 
d 2 a/ir dp\ dyi ab for p(770)°, i^*(892)° and 0(1020) mesons in the rapidity range \yi ab \ < 1.0 in 
bins of pt- The first error is statistical and the second systematic. For each bin in pr the range 
as well as the bin-centred value p^ n are given. 





dcr/dyi ab [fib] 


yiab 


P° 


(K*° + K*°)/2 


<P 


[-1.0,-0.5] 


11.0 ±1.0 ±1.2 


3.36 ±0.35 ±0.72 


1.44 ±0.25 ±0.22 


[-0.5,0.0] 


13.1 ± 1.1 ± 1.4 


2.52 ±0.27 ±0.36 


1.08 ±0.12 ±0.16 


[0.0,0.5] 


10.4 ± 1.5 ± 1.1 


3.07 ±0.30 ±0.44 


1.44 ±0.13 ±0.22 


[0.5, 1.0] 


14.6 ±1.3 ±1.5 


4.28 ±0.44 ±0.79 


1.61 ±0.33 ±0.25 



Table 2: Inclusive non-diffractive photoproduction single differential cross sections da jdyi^ 
for p(770)°, i^*(892)° and 0(1020) mesons in the transverse momentum range pr > 0.5 GeV 
in bins of yi a b- The first error is statistical and the second systematic. 
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(K*° + K*°)/2 




7 p(Hl) 


(da/dy lab ) lyiabl<1 [/xb] 


23.6 ±2.7 


5.22 ±0.60 


1.85 ±0.23 


T [GeV] 


0.151 ±0.011 


0.166 ±0.012 


0.170 ±0.012 


(£ T ) [GeV] 


1.062 ±0.018 


1.205 ±0.020 


1.333 ±0.022 


(£^ n ) [GeV] 


0.287 ±0.018 


0.313 ±0.020 


0.314 ±0.022 


(pr) [GeV] 


0.726 ±0.027 


0.810 ±0.030 


0.860 ±0.035 


pp (STAR) 


(p T ) pp [GeV] 


0.616 ±0.062 


0.81 ±0.14 


0.82 ±0.03 


Au-Au (STAR) 


MauAu [GeV] 


0.83 ±0.10 


1.08 ±0.14 


0.97 ±0.02 



Table 3: The parameters (da / dy iab ) and T = E To /n for p°, K*° and mesons from 

a fit of function (1) to the differential cross sections. The average transverse energy (Et), 
kinetic energy (Elp n ) and momentum (p T ) are also presented. The errors correspond to the 
quadratically summed statistical and systematic errors. Also shown are measurements in pp 
and Au-Au interactions at nucleon-nucleon centre-of-mass energy ^/snn = 200 GeV [4] at 
central rapidities. 



Experiment 


Measurement 


R(<f)/K* ) 


HI 


7 p, (W) = 210 GeV, M < 1 


0.354 ±0.060 


STAR 


pp, v/s = 200 GeV, \y\ < 0.5 


0.35 ±0.05 


Au-Au, ^/s NN = 200 GeV, \y\ < 0.5 


0.63 ±0.15 



Table 4: The ratio R(4>/K*°) of the total cross-sections for (f> and A'* production obtained 
in 7p collisions (HI) at (W) = 210 GeV. The errors correspond to the quadratically summed 
statistical and systematic errors. Also shown are measurements in pp and Au-Au interactions at 
nucleon-nucleon centre-of-mass energy y/s NN = 200 GeV [4] at central rapidities. 
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Figure 1: The unlike-sign di-pion mass spectrum after subtracting the like-sign contribution, 
normalised to the total number of entries. The solid and dashed curves correspond to the 
PYTHIA simulation with and without Bose-Einstein correlations (BEC), respectively. 
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Figure 2: The invariant mass spectra for ix + ix~ in a) and 6), for K ± tt t in c) and for K + K~ in ci). 
The full curves show the result of the fit; the dashed curves correspond to the contribution of the 
combinatorial background B(m). In b), the data and the fit F{m) are shown after subtraction of 
the combinatorial background B(m); the dotted and dash-dotted curves show the contributions 
from uj and K* reflections, respectively. In c), the dotted curve corresponds to the contribution 
of the reflections and the dash-dotted curve corresponds to the contribution of the A'* (892) 
signal. In d), the dotted curve corresponds to the contribution of the 0(1020) signal. 
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Figure 3: The inclusive differential non-diffractive cross sections for p°(770), J^*°(892) and 
0(1020) mesons measured in a) as a function of transverse momentum for \yi a b\ < 1 and in b) 
as a function of rapidity for p T > 0.5 GeV. The curves on the figure a) correspond to the power 
law, equation (1), with n = 6.7. The ratios of data to Monte Carlo predictions "Data/MC" are 
shown for the PYTHIA (full points) and PHOJET (empty points) simulations as a function of 
transverse momentum for \yi a b\ < 1 in c) and as a function of rapidity for pr > 0.5 GeV in d). 
Statistical and systematic errors are added in quadrature. 
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Figure 4: The inclusive invariant differential cross sections as a function of transverse momen- 
tum. The curves show the results of fits to the power law, equation (1). Statistical and systematic 
errors are added in quadrature. 
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